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The spin dynamics of the S = 1 Ni-chains in mixed-spin antiferromagnets P^BaNiOs and Nd I Y2-iBaNiOs is described in 
terms of a simple Ginzburg-Landau Lagrangian coupled to the sublattice of rare-earth ions. Within this framework we obtain a 
theoretical explanation for the experimentally observed coexistence of Haldane gap excitations and long-range magnetic order, 
as well as for the increase of the Haldane gap energy below the Neel point. We also predict that the degeneracy of the Haldane 
triplet is lifted in the magnetically ordered phase. The growth of the transverse gap is shown to follow from the magnon 
repulsion. The theoretical results are consistent with the available experimental data. 



Much excitement in the field of low-dimensional quan- 
tum magnetism was caused by the theoretical work of 
Haldane, who predicted an energy gap in the excitation 
spectrum of a 1-dimensional Heisenberg antiferromagnet 
(1-D HAF) with integer spin jij. At the beginning this 
came as a surprise, since for the half-integer-spin 1-D 
HAF the excitation spectrum was known to be gapless. 
Very soon Haldane's conjecture was confirmed by many 
numerical and experimental studies Q , and the "Haldane 
gap" phenomenon is by now rather well understood. A 
very challenging outstanding problem now is to study the 
quasi- 1-D case, when 3-D magnetic interactions are suffi- 
ciently strong to destroy the Haldane singlet and produce 
long-range magnetic order at low temperature, yet some 
purely quantum-mechanical effects are preserved, thanks 
to the dominance of 1-D interactions. Much experimen- 
tal and theoretical work in this direction was done on 
CsNiCl3 || and isostructural integer-spin compounds [Q . 
The result was a profound understanding of the connec- 
tion that exists between Haldane excitations and "nor- 
mal" spin waves. In CsNiCLj for example, the triplet of 
1-D Haldane gap modes is readily observed above the 
Neel temperature Tjy. As TV is approached from above, 
the gap vanishes at the 3-D magnetic zone-center, driv- 
ing a soft-mode transition. In the ordered phase two of 
the three modes become conventional gapless spin waves, 
while the third "longitudinal" mode, not accounted for 
by the linear spin wave theory, persists as a gapped ex- 
citation. 

Recently in a series of inelastic neutron scattering stud- 
ies of quasi- ID mixed-spin antiferromagnets with the gen- 
eral formula i? 2 BaNi0 5 {R =Pr, Nd or Nd x Yi_ x ) Zhe- 
ludev et al. H demonstrated that the scenario realized 
in CsNiCi3 is by no means universal. In Ji^BaNiOs com- 
pounds 3-D magnetic ordering produces finite static mo- 
ments on both the 5=1 Ni 2+ and the i? 3+ magnetic 
ions H , yet the 1-D gap excitations propagating on the 
Ni-chains are only weakly affected: they persist in the 
ordered phase as gapped modes and coexist with conven- 



tional acoustic spin waves all the way down to T = 0. 
These excitations show no signs of softening at the tran- 
sition point at any wave vector. In the entire phase di- 
agram they have a purely 1-D dynamic structure factor 
that is practically indistinguishable from that of Haldane 
excitations in Y2BaNiOs [Q, a "clean" Haldane-gap ma- 
terial, where only the Ni-sites are magnetic, and no long- 
range ordering occurs. In ii^BaNiOs the only significant 
(and rather unexpected) effect of 3-D ordering is that 
the gap energy increases in the ordered phase, the in- 
crease being roughly linear with (Tj\r — T) . The survival 
of 1-D quantum spin excitations in the ordered state is 
indeed remarkable since the transition temperatures are 
rather high, and are quite comparable to the gap en- 
ergy. The purpose of this letter is to explain the coexis- 
tence of long-range order and Haldane gap excitations in 
i?2BaNi05 materials. We propose that, contrary to the 
case of CsNiCLj, in these compounds Ni-chains do not in- 
teract between themselves, but instead each is coupled to 
the sublattice of magnetic i?-ions. We derive the exper- 
imentally observed growth of the gap below TV using a 
simple Ginzburg-Landau Lagrangian J8| . In addition, we 
predict that the degeneracy of the Haldane excitations is 
partially removed in the ordered phase and notice an in- 
triguing analogy between this gap splitting in i?2BaM05 
materials and the famous Higgs mechanism of lifting the 
mass degeneracy in particle physics. 

Integer spin chains are traditionally described in terms 
of the Lagrangian of the Non-Linear Sigma Model 
(NLSM). In this approach one approximates the field of 
spin site-operators Si (Sf — s(s + 1)), retaining only the 
Fourier components near k x = 0, and n Q. This is ac- 
complished by a change in real space variables 

£ = S (-i)V«+r(i) , (i) 

where both fields <p(i), and change slowly on the scale 
of one lattice spacing. According to this definition, (p{i) 
is the unit vector in the direction of local staggered mag- 
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netization, and l(i) is the component of local magne- 
tization perpendicular to <p(i) (<p(i) ■ l(i) — 0). In the 
limit of large s three components of cp commute with each 
other, while the commutation rules of I remain nontriv- 
ial. It can be shown that the correct dynamics for the 
new variables follows from a NLSM Lagrangian: 
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where ip is subject to the constraint ip 2 — 1. Here 
g = 2/s is a dimensionless coupling constant, measuring 
the strength of quantum fluctuations, and v = 2Js is the 
spin wave velocity. The magnetic moment / is given by 
I = 1 /gv ip x dip/dt. This mapping of the 1-D HAF to the 
NLSM (||) becomes exact for integer s — ► oo, but it gives a 
meaningful and qualitatively correct approximation even 
for s=l. From this description it can be derived that 
in 1+1 dimensions any strength of quantum fluctuations 
g is sufficient to destroy long range correlations between 
spins. A finite correlation length £ is always accompa- 
nied by a gap A = v/£ in the spin excitation spectrum. 
For half-integer 1-D HAF the topological term not in- 
cluded in the Lagrangian (Q) prevents the appearance of 
the finite correlation length, and spin correlations remain 
scale-free. But for integer s this gap, commonly known 
as the Haldane gap, is present. It has a small g (large 
s) asymptotics of A ~ exp(— 2it/g) = exp(— its). For the 
s = 1 HAF the Haldane gap was found numerically to be 
A = 0.41J g. 

The non-linear sigma model is not very convenient 
for practical calculations. To simplify things further we 
follow the approach of Affleck [||, and replace the La- 
grangian of NLSM with that of the quantum Ginzburg- 
Landau model 
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where for convenience we have simultaneously changed 
variables to <p = 1/y/g = y/s/2 (p. The new La- 
grangian follows from (0) if the constraint p 2 = 1 is 
relaxed, making the vector of staggered magnetization 
soft, and phenomenological quadratic and quartic terms 
are introduced. The coefficient in front of the quadratic 
term is selected to reproduce the correct value of the gap 
(in NLSM this gap is generated dynamically), while the 
only condition imposed on the quartic term is that A > 0, 
ensuring the overall stability for large \<p\. The physical 
meaning of Eq. (JsJ) is rather transparent. It describes the 
propagation of a triplet of magnons (excitations of local 
staggered magnetization) with a given spin wave velocity 
v, and gap A. The quartic term describes the repulsive 
magnon-magnon interaction. 



The Lagrangian (|3|) was successfully used to qualita- 
tively explain the behavior of coupled ID spin chains in 
CsNiCl3 [0j . In this case the i-th chain is described by Eq. 
(0) with <p <pi, while the interchain coupling is given 
by —2s J' '(pi ■ (pi+\. On the mean-field level the behav- 
ior of the excitations is the following: the gap is reduced 
by the interchain coupling to A 2 ^ = A 2 (T) — AvZsJ', 
where Z is the number of chains coupled to a given 
chain. The 0(3) symmetry is spontaneously broken for 
A 2 (T) < AZsvJ'. In this case \(p\ has a nonzero ex- 
pectation value of y/(4vZs,J' - A 2 (T))/4uA. The gap 
naturally grows with temperature and eventually wins 
over the interchain exchange, restoring the 0(3) symme- 
try. The Neel temperature is determined by the condi- 
tion A 2 (Tjy) = AvZsJ'. It is important to understand 
what happens to the triplet of magnons below Tjy. As 
Tn is approached from above, the gap in all three exci- 
tations at the positions of future magnetic Bragg peaks 
decreases and hits zero at Tn- Below Tjy two excita- 
tions corresponding to the two Goldstone modes of the 
Neel order parameter remain gapless. At the same time 
the longitudinal branch re-acquires the gap along with a 
finite lifetime due to decay into two gapless modes. 

The "traditional" scenario || described above is in 
agreement with experimental findings in CsNiCl3. How- 
ever, for i?2BaNi05 we can expect a different picture. 
The structure of these materials is such that there are 
no strong direct superexchange paths between the Ni 
chains, which explains why in isostructural Y2BaNi05 
no magnetic ordering is observed 0. On the other 
hand, magnetic rare earth ions substituted for Y are cou- 
pled between themselves. To a good approximation one 
can assume that these sites form a separate sublattice 
which orders at rather high temperatures (Tn = 48K for 
Nd2BaNi05). The coupling J' between the Ni and R sub- 
lattices is weak, yet finite. The finite ordered moments 
on the Ni-sites may be interpreted as the singlet-ground- 
state Ni-chains being polarized by a staggered exchange 
field from the magnetically ordered sublattice of i?-ions. 
To describe this effect an "external field" term should 
be added to the Ginzburg-Landau Lagrangian (^). It is 
given by J'n-\[2s(p = H s ■ <p, where ft is the staggered mo- 
ment of the R sublattice. Such a term explicitly breaks 
the 0(3) symmetry and causes (p to acquire a nonzero 
expectation value along H s . As a consequence, the de- 
generacy of the triplet of Haldane gap excitations is par- 
tially lifted. The energy gap in the magnon, polarized 
parallel to H s , is different from that in the two transver- 
sal branches. Both gaps can be derived from V - the 
potential energy density of the Lagrangian: 
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The usual formula for this is A 2 /w = d 2 V /d(p 2 a . It gives 
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A^=Al + l2vA<f>% 
A 2 ± = A„ + 4uA0g 



(5) 



Here fa is the expectation value of \<p\ determined by 
minimization of (E) through 



H s = A 2 /v cf> + 4A0;j 
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or, approximately, 4>o — vH s /A^. Note that both energy 
gaps are increased compared to their value in the absence 
of the staggered field. This observation relies on the fact 
that A > 0, i.e. a repulsive magnon-magnon interaction. 
It is precisely the positive A, which guaranties the sta- 
bility of Eq. (^J) and of the Lagrangian (^) itself, and, 
therefore, should be satisfied. 

As we mentioned in the introduction, the Lagrangian 
of NLSM is traditionally used to describe the behavior of 
quantum spin chains. The effect of the staggered mag- 
netic field in NLSM , coupled to the staggered magnetiza- 
tion via the term ft ■ h/g, was studied by Nelson and Pel- 
covits [ fL0| using renormalization group methods. Their 
RG flow indicates that indeed the gap increases (the cor- 
relation length decreases) in strong enough fields. The 
one loop corrections calculated by them do not give a def- 
inite answer about the behavior of the correlation length 
in weak fields. Using our results for the Ginzburg-Landau 
Lagrangian we can conclude that the correlation length 
decreases in this case as well. 

Our predictions for the energy gap can be qualita- 
tively understood by working with the original Heisen- 
berg Hamiltonian and treating the staggered field as a 
weak perturbation. Let us denote the Haldane ground 
state by \G) and label the lowest-energy triplet excita- 
tions at q = 7r as |J5,+1), \E,0) and \E,—1), accord- 
ing to their value of S z . Let us now consider the rel- 
evant matrix elements of the staggered field operator 
H s = H s J2j(~^ySj- These are non-zero only between 
states with momentums differing by exactly 7r. In par- 
ticular, (G\H S \G), as well as the matrix elements be- 
tween any two of the \E, +1), \E,0) or \E, — 1) states 
are strictly zero. As a result, there are no first-order 
corrections to the spin gap at q = ir. To calculate the 
2-nd order correction 5Eq to the ground state energy 
we can use the well-known fact that the triplet of Hal- 
dane modes pretty much exhausts the spectral weight 
at q = ir. Thus, only the mixing of \G) with single- 
particle excitations needs to be considered. Since H s 
conserves the z-component of total spin, its only non- 
zero matrix element is the one between \G) and \E, 0). 
The correction to the ground state energy is thus given 
by 6E G = _ K^ |g°l G >l 2 ee -e. We now proceed to calcu- 
late the energy corrections for the three excited states at 
q = ir. Two effects need to be considered: i) their mixing 
with the ground state and ii) their mixing with 2-particle 
excited states. Let us first assume that there is no repul- 
sion between magnons, and that doubly-excited states 



at q = 0, which we label as \2E,+1,0), \2E,—l, 0), and 
\2E,0,0), are produced by adding two non-interacting 
q = ir excitations to the system. Second order corrections 
to the energies of singly-excited states are then given by: 
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correction to the gap energies then become zero for the 
transverse gap Aj_, and +e oc H 2 for the longitudinal 
gap Am. Magnon repulsion effectively reduces the nega- 
tive second-order corrections to the single-particle excita- 
tion energies that are due to mixing with double-excited 
states. Even for the transverse modes magnon repulsion 
will thus lead to an increase of the gap energy. Both 
Aj_ and A|| will increase quadratically with iT s , though 
for A|| the effect will be more pronounced. This is to- 
tally consistent with results derived above, where the 
positive parameter A represents repulsive interactions be- 
tween magnons. 

An increase of the energy gap below the transition 
point is exactly what was observed in neutron scatter- 
ing experiments on Nd2BaNiOs, Nd x Y2_ x BaNi05, and 
Pr 2 BaNi0 5 §. The data were taken on powder samples 
or single crystals of small size, which obviously made re- 
solving the mode splitting below TV impossible. In the 
powder measurements at least, the observed gap always 
corresponds to the lowest-energy mode. In Fig. 1 we 
show the combined data from all experiments, plotting 
the increase of the square of gap energy (relative to that 
measured in Y2BaNi05 |ll]] at the same temperature) 
against the square of staggered magnetization of the Ni 
sublattice. The data are consistent with the predicted 
linear dependence. Moreover, data points for different 
compounds seem to collapse onto a single curve, even 
though the Neel temperatures vary considerably. The 
mechanism is therefore insensitive to the details of ex- 
change in the i?-sublattice, and between Ni and R sub- 
lattices. 

In a much obscured form the effect of staggered field 
has been previously observed in the well-known Haldane- 
gap material NENP [Q. In this compound applying 
an external uniform field H produces a weak effective 
staggered component due to some special structural fea- 
tures. On top of a very pronounced linear splitting of 
the Haldane triplet by the uniform field for each mode 
one observes a slight energy shift that is positive and 
quadratic with H . This shift was attributed to the weak 
staggered field component that for NENP is insepara- 
ble from the uniform one |l3f| . Numerical calculations 
have confirmed this conclusion |l4|| . Note that unlike in 
NENP, in i?2BaNi05 systems we are able to see the effect 
of staggered field much more clearly, since the dominant 
uniform field component is absent. 

It is interesting to mention that the effect of gap in- 
crease and splitting due to interaction with another sub- 
lattice has a close analogue in the field of particle physics. 
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It is very similar to the Higgs mechanism by which the 
mass degeneracy in a multiplet of elementary particles 
is lifted. As in our case, the splitting of masses (en- 
ergy gaps relative to a vacuum of particles) is a result of 
spontaneous symmetry breaking, produced by the inter- 
action with Higgs particles. The three Ni chain Haldane 
excitations of different polarizations play the role of a 
multiplet of particles of initially equal masses (quarks for 
example), while the conventional acoustic spin waves, re- 
ferred to as mixed i?-Ni excitations by Zhcludcv et al, 
play the role of the elusive Higgs particle, which is mass- 
less in our case. This analogy should not be taken too 
seriously due to the differences in Hamiltonians: it is the 
SU(2) symmetry which is spontaneously broken by the 
Higgs mechanism, and the term describing the coupling 
between matter and Higgs particles is more complicated 
than a simple —H s <p term in our case. Nevertheless, the 
essence of the effect is the same: spontaneous symmetry 
breaking caused by interaction with some external field 
lifts the degeneracy of masses (gaps) . 

In conclusion, we have provided a simple theoretical 
explanation for the coexistence of long range magnetic 
order and Haldane gap excitations in i^BaNiOg systems. 
The Ni-chain gap modes are described in terms of a sim- 
ple Ginzburg-Landau Lagrangian. The increase of the 
gap below the Neel temperature, previously observed ex- 
perimentally in several compounds, immediately follows 
from this approach. The qualitatively new prediction is 
that the degeneracy of the Haldane triplet is lifted be- 
low the transition point. The magnon polarized along 
the vector of staggered magnetization has a larger gap 
than the two others magnon branches. Experimental 
tests of this prediction, including inelastic neutron scat- 
tering experiments on large high-quality single samples, 
are planned for the near future. 

We thank Victor Emery and Igor Zaliznyak for help- 
ful discussions, and V.N. Muthukumar for bringing Ref. 

to our attention. This work was carried out under 
Contract No. DE-AC02-76CH00016, Division of Material 
Science, U.S. Department of Energy. 



Phys. Condens. Matter 3, 715 (1991). 

[4] A. Harrison, M. F. Collins, J. Abu-Dayyeh, and 
C. V. Stager, Phys. Rev. B 43, 679 (1991); M. Enderle 
et al, Europhys. Lett. 25, 717 (1994). 

[5] A. Zheludev, J. M. Tranquada, T. Vogt, and D. J. Buttrey, 
Europhys. Lett. 35, 385 (1996); Phys. Rev. B 54, 6437, 
and 7210 (1996); T. Yokoo, A. Zheludev, M. Nakamura, 
and J. Akimitsu, Phys. Rev. B 55, 11516 (1997). 

[6] D. J. Buttrey, J. D. Sullivan, and A. L. Rheingold, J. Solid 
State Chem. 88, 291 (1990); V. Sachan, D. J. Buttrey, 
J. M. Tranquada, and G. Shirane, Phys. Rev. B 49, 9658 
(1994); A. Zheludev, J. P. Hill and D. J. Buttrey, Phys. 
Rev. B. 54, 7216 (1996). 

[7] J. Darriet and L. P. Regnault, Solid State Commun. 86, 
409 (1993); J. F. DiTusa et al, Physica B 194-196, 181 
(1994); G. Xu et al, Phys. Rev. B 54, R6827 (1996). 

[8] I. Affleck, Phys. Rev. Lett. 62, 474 (1989); I. Affleck, and 
G. F. Wellman, Phys. Rev. B 46, 8934 (1992). 

[9] O. Golinelli, T. Jolicoeur, and R. Lacaze, Phys. Rev. B 45, 
9798 (1992). 

[10] D. R. Nelson and R. A. Pelcovits, Phys. Rev. B 16, 2191 
(1977). 

[11] T. Sakaguchi, K. Kakurai, T. Yokoo, and J. Akimitsu, J. 

Phys. Soc. Jpn. 65, 3025 (1996). 
[12] M. Chiba et al., Phys. Rev. B 44, 2838 (1991). 
[13] T. Sakai and H. Shiba, Jrn. Phys. Soc. Jpn. 63, 867 (1994). 
[14] P. P. Mitra and B. I. Halperin, Phys. Rev. Lett. 72, 912 

(1994). 

FIG. 1. Experimentally measured increase of the square 
of the gap energy in .ffeBaNiOs samples relative to that 
in Y2BaNi05 at the same temperature, plotted against the 
square of staggered magnetization of Ni sites. 



[1] F. D. M. Haldane, Phys. Lett. 93A, 464 (1983); Phys. Rev. 
Lett. 50, 1153 (1983). 

[2] For a comprehensive list of references see some recent pub- 
lications, e. g., L. P. Regnault, I. Zaliznyak, J. P. Renard, 
and C. Vettier, Phys. Rev. B 50, 9174 (1994); A. Zheludev 
et al, Phys. Rev. B 53, 15004 (1996); G. Xu et al, Phys. 
Rev. B 54, R6827 (1996). 

[3] See for example W. J. L. Buyers et al, Phys. Rev. Lett. 
56, 371 (1986); R. M. Morra, W. J. L. Buyers, R. L. Arm- 
strong, and K. Hirakawa, Phys. Rev. B 38, 543 (1988); 
K. Kakurai, M. Steiner, R. Pynn, and J. K. Kjems, J. 



4 




( 2 A9W) V Z V 



